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Abstract 

We use the coupled cluster method for infinite chains comple- 
mented by exact diagonalization of finite periodic chains to discuss 
the influence of a third-neighbor exchange J3 on the ground state of 
the spin-^ Heisenberg chain with ferromagnetic nearest-neighbor in- 
teraction J\ and frustrating antiferromagnetic next-nearest-neighbor 
interaction J2. A third- neighbor exchange J3 might be relevant to 
describe the magnetic properties of the quasi-one-dimensional edge- 
shared cuprates, such as L1VCUO4 or LiCu202- In particular, we cal- 
culate the critical point J| as a function of J3, where the ferromagnetic 
ground state gives way for a ground state with incommensurate spi- 
ral correlations. For antiferromagnetic J3 the ferro-spiral transition 
is always continuous and the critical values J| of the classical and 
the quantum model coincide. On the other hand, for ferromagnetic 
J3 < —(0.01 . . . 0.02) I Ji I the critical value J% of the quantum model 
is smaller than that of the classical model. Moreover, the transition 
becomes discontinuous, i.e. the model exhibits a quantum tricritical 
point. We also calculate the height of the jump of the spiral pitch 
angle at the discontinuous ferro-spiral transition. 
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PACS codes: 

75.10.Jm Quantized spin models 

75.45.+j Macroscopic quantum phenomena in magnetic systems 
75.10.-b General theory and models of magnetic ordering 



1 Introduction 

The recent observation of spiral (helical) magnetic ground states in sev- 
eral chain cuprates, such as LiVCuO.4, LICU2O2, NaCu202, Li 2 ZrCu04, and 
Li 2 Cu02 [IHTT]. which were identified as quasi-one- dimensional (ID) frus- 
trated spin- 1/2 magnets with ferromagnetic (FM) nearest-neighbor (NN) 
in-chain J\ and antiferromagnetic (AFM) next-nearest-neighbor (NNN) in- 
chain interactions J2 has stimulated intensive investigations of frustrated ID 
Heisenberg ferromagnets, see, e.g., Refs. [T2~H2"T] . The ID J\- J% model consid- 
ered in the most theoretical papers may serve only as the minimal model to 
describe the magnetic properties of these materials. Several extensions, such 
as exchange anisotropy [HJI2I] or interchain coupling [T2JI22H22] might be 
relevant to explain experiments. In addition to the NN and NNN exchange 
integrals, J\ and J2, also an exchange coupling to 3rd neighbors, i.e. J3, or 
even couplings to farther distant neighbors could play a role in real materials. 
One mechanism to induce such exchange interactions is strong spin-phonon 
interaction for frustrated chains within the antiadiabatic limit [2B]. Even 
if these additional couplings are small, their influence on the spiral ground 
state (GS) correlations might be noticeable. In particular, there is a signifi- 
cant influence of J3 on the critical frustration J| at which the FM GS gives 
way for the spiral GS, see below. Except the possible relevance of the J\-J<r 
J3 model for real materials the consideration of such a model is interesting 
in its own right as a basic model to study frustration effects in ID quantum 
spin systems. Moreover, very recently it has been argued that the magnetic 
properties of the kagome-like mineral volborthite Cu3V207(OH)2-2H20 can 
be described by an effective chain model with farther distant frustrating ex- 
change couplings [27]. The corresponding general Heisenberg Hamiltonian 
H with NN exchange Ji, NNN exchange J 2 and farther in-chain exchange 
interactions J n reads 

+ ... • (1) 
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Motivated by the experiments on the edge-sharing chain cuprates we focus 
on the spin-1/2 J1-J2-J3 model with FM Ji and frustrating AFM J 2 > 0. To 
the best of our knowledge this model has been investigated so far only in an 
early paper of Pimpinelli et al. [12] using spin-wave theory. 

Here we use the coupled cluster method (CCM) for infinite chains com- 
plemented by exact diagonalization (ED) of finite chains (periodic bound- 
ary conditions imposed) to investigate spiral GS correlations. Both meth- 
ods have been successfully applied to study the spiral ordering of the J1-J2 
model [inmHEH] . In Refs. [131115] it was demonstrated that the CCM results 
are in good agreement with the DMRG data. However, in order to take 
into account the J 3 bonds properly, we go beyond the so-called SUB2-3 ap- 
proximation used in Refs. [131112] and consider an improved approximation, 
namely the LSUB4 approximation, see below. 

The paper is organized as follows. In Sec. [2] discuss briefly the classical 
GS. In Sec. |3j we provide a brief illustration of the CCM and describe its 
application on the considered model. In Sec. [Hour results for the FM-spiral 
phase transition and for the pitch angle in the spiral phase are presented and 
discussed. In Sec. H2 we summarize our findings. 

2 The classical model 

First we discuss the GS of the classical model (spin quantum number s — > 00). 
For the usual Ji-J 2 model, i.e. the model with J n = 0, n > 3, studied in 
many papers the critical frustration J 2 is Jf = |Ji|/4. For J 2 > J| there is a 
spiral GS with a canting angle between NN (pitch angle) 7 given by cos 7 = 
|Ji|/(4J 2 ). This helix interpolates between a FM chain at < J 2 < J 2 and 
two decoupled AFM chains at J 2 /| J\\ = 00. Noteworthy, J| is unaffected by 
quantum effects, see, e.g., Refs. [T2 | [T3 | IT5 | [T9]. 

If farther AFM couplings J n (n > 3) are relevant the destabilzation of 
the FM GS sets in for smaller values of J 2 . Extending the classical model 
including arbitrary J n 7^ (n > 3) one finds for the critical NN exchange 



This expression has been derived assuming a continuous transition between 
the FM and the spiral GS's. It holds for arbitrary AFM long-range couplings 
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J n with n > 3. In case that some of the exchange couplings are FM, i.e. 
J n < for certain n > 3, the Eq. (J2]) holds only if the AFM couplings 
dominate. Assuming J 2 > we find as the criterion for the validity of 
Eq. © 



J 2 > 



1 00 



12 

n=3 



or equivalently 



w>4E na ( na - 4 ) J - ( 4 ) 

n=3 

If this condition is violated, in crossing the FM-spiral phase boundary, the 
spiral GS "jumps" from a finite pitch angle 7^ 7^ to 7 = in the FM GS. 
For the simplest case J 3 7^ and J n = (n > 3) the classical model was 
considered by Pimpinelli et al. [12]. They found for the critical frustration 
J 2 in case of continuous transition J| = (|J;iJ — 9J 3 )/4 in accordance with the 
general expression ([2]). The classical pitch angle in the spiral phase is given 

by cos 7 = ^— J 2 + \fj\ + 3J 3 (3J 3 + \Ji\)j /6J3. According to the general 
Eqs. ([3]) and (j4j) for the J1-J2-J3 model the classical FM-spiral transition is 
discontinuous at J 3 < — ^ or equivalently at J 3 < — y. 

For J 3 < and J n = 0, n > 3, the height of the jump of the pitch angle 
7t at the transition is given by [12] 



cos 7 t = - Ji ( 2 J 3 +V4J| + 4 JsJ,) -- (5) 
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or equivalently 



008^ = — --. (6) 

Considering other simplified classical models with a single FM long-range 
coupling J no , i.e. a model with J\ < 0, J 2 > 0, J no < 0, J n = {n > 2 
and ri 7^ no), the classical discontinuous FM-spiral transitions occurs ac- 
cording to Eqs. (j3J) and (T4]) at — J no > 12J 2 / (uq (uq — 1)) (or equivalently 
at — J no > 3 1 A|/(^o( n o — 4)). In other words, even a tiny but fairly long 
range ferromagnetic coupling may introduce a discontinuous behavior at the 
critical point. 
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3 The Coupled Cluster Method (CCM) 



For the sake of brevity, we will outline only some important features of the 
CCM which are relevant for the model under consideration. The interested 
reader can find more details concerning the application of the CCM on the 
frustrated Heisenberg magnets with non-collinear GS's in Refs. [T3 | IT9" | I29T 
[37]. For more general aspects of the methodology of the CCM, see, e.g., 
Refs. HEM- 

First we mention that the CCM approach yields results in the thermody- 
namic limit N — > oo. The starting point for a CCM calculation is the choice 
of a normalized reference (or model) state |<3>). Related to this reference 
state we then define a set of mutually commuting multispin creation oper- 
ators Cj~, which are themselves defined over a complete set of many-body 
configurations I. For the considered frustrated spin system we choose a spiral 
reference state with spiral spin orientations along the chains (i.e., pictorially, 
1^) = I " ' ' )) characterized by a pitch angle 7, i.e. |$) = |$(7)). 

Such states include the FM state (7 = 0) as well as the Neel state (7 = 7r). 
Next, we perform a rotation of the local axis of the spins such that all spins 
in the reference state align along the negative z axis. This rotation by an 
appropriate local angle 5 n of the spin on lattice site i is equivalent to the 
spin-operator transformation 

s x n = cos 5 n s x n + sin 5 n s z n ; s y n = s v n ; s z n = - sin 5 n s x n + cos 5 n s z n , (7) 

where s^, s^, are the spin operators in the rotated coordinate frame. The 
local rotation angle 5 n is related to the pitch angle 7 of the spiral reference 
state by 5 n = 717. In this new set of local spin coordinates the reference state 
and the corresponding multispin creation operators Cf are given by 

|$) = I 444 • ■ ■ ); Cf = sf , sfs^ , sfs^sf , . . . , (8) 

where the indices n,m,k, . . . denote arbitrary lattice sites. In the rotated 
coordinate frame the Hamiltonian becomes dependent on the pitch angle 7. 
It reads 

m=l n 

+ [cos(m 7 ) - 1] (s+s+ +m + KK+m) + 2 sin(m 7 ) [s+s z n+m 

~ s n § n+m + s n s n+m — s n s n+m ] + 4 cos(m7)s n s n+m , (9) 
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where sj n = s x i n ± is\ n . 

With the set |J$),C/} the CCM parametrization of the exact ket GS 
eigenvector \^f) of the many-body system is given by 

|M/)=e 5 |$>, S = J2 a iCf- (10) 

The CCM correlation operator S contains the correlation coefficients aj, 
which can be determined by the so-called set of the CCM ket-state equations 

(<S>\Cye- s He s \<!>} =0 ; VI ^ 0, (11) 

where Cj = (Cj~) . Each ket-state equation belongs to a specific creation op- 
erator Of = s+, s+s+ , s+s+s^ , • • • , i.e. it corresponds to a specific set (con- 
figuration) of lattice sites n,m,k, . . . . By using the Schrodinger equation, 
H\^J) = E\ty), we can write the GS energy as E = (&\e~ s He s \&) = E(j), 
which depends (in a certain CCM approximation, see below) on the pitch 
angle 7. 

In the quantum model the pitch angle may be different from the correspond- 
ing classical value 7 c i. Therefore, we do not choose the classical result for the 
pitch angle in the quantum model, rather, we consider 7 as a free parameter 
in the CCM calculation, which has to be determined by minimization of the 
CCM GS energy £(7), i.e. dE/d-f\ 1=lqu = 0. 

For the many-body quantum system under consideration it is necessary to use 
approximation schemes in order to truncate the expansions of S in Eq. f llOp 
in a practical calculation. In Refs. [13] and [TH] it has been demonstrated, 
that for the J1-J2 model the so-called SUB2-3 approximation leads to results 
of comparable accuracy to those obtained using the DMRG method [H]. In 
this approximation all configurations are included which span a range of no 
more than 3 contiguous sites and contain only 2 or fewer spins. Taking into 
account the J3 bond we have to extent this approximation in order to take 
into account configurations including a range of 4 contiguous sites. The cor- 
responding approximation is the so-called LSUB4 approximation, see e.g., 
Refs. [35l|39lll0] . Within this approximation multispin creation operators of 
one, two, three or four spins distributed on clusters of four contiguous lattice 
sites are included. 

In addition, for the determination of the quantum tricitical point we 
have also used higher LSUBn approximations, see Sec. HI However, the 
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numerical complexity increases tremendously, since (i) the number of ket- 
state equations ( II ip increases exponentially with n, (ii) there are two free 
parameters J2, J3 which have to be varied by very small increments to find 
the transition points, and (iii) the determination of the quantum pitch angle 
7 gu itself requires the iterative minimization of the ground state energy for 
each set of J 2 , J3. Hence, except for the determination of the quantum 
tricitical point, we have restricted our CCM calculations to the CCM-LSUB4 
approximation. 



4 Results for the spin- 1/2 quantum model 

In what follows we set Ji = — 1 if not stated otherwise explicitly. The 
phase diagram of the J1-J2-J3 chain with FM J\ = —1 obtained by the 
CCM and by the ED is shown in Fig. [TJ For the classical as well as for the 
quantum model the transition from the FM to the spiral GS can be second 
(J 3 > J 3 ) or first order (J 3 < jj), i.e., the pitch angle 7 does (first order) 
or does not (second order) jump from 7 = to a finite value 7 = 77- > 
at the transition. The tricritical point T = {J 2 ,J 3 ), i-e. that point at 
the transition line where the second-order transition goes over in a first- 
order transition is T c \ = (jJ oi , J 3 cl ) = (§,— -h) for the classical model [T2] . 
see the black square in Fig. [TJ Due to quantum fluctuations this point is 
shifted to smaller values of | J 3 \. The LSUB4-CCM estimate for the quantum 

.„,,„, rpLSUBi T LSUBi 

tricritical point is T^ um = {J 2 qu , J 3 qu ) = (0.283,-0.013), see the 
red square in Fig. [TJ This result is quite close to the spin-wave result [12] 

Tq™ — (^2 9 " 1 ') = (0.25,0). We have determined the quantum tricrit- 
ical point also using higher CCM-LSUBn approximations, namely LSUB6, 

rpLSUBS rpLSUBd 

LSUB8 and LSUB10. The corresponding results are {J 2 qu iJz qu ) — 

rpLSUBS, rpLSUBS 

(0.274,-0.011), (J 2 qu ,J 3 qu ) 

rpLSUBW rpLSUBlO 

= (0.268,-0.009), and {J 2 qu , J 3 qu ) = (0.266,-0.007). Obviously, 
these values are quite close to each other, and there is a slight shift towards 
{J 2 qu , J 3 qu ) of Ref. [TJ when increasing the order of CCM approximation. 

The ED data for the transition line are in very good agreement with the 
CCM data. There is only a slight deviation for stronger FM J 3 bonds, i.e. 
f° r ^3 ^ —0.2. Note, however, that in a tiny interval —0.1 < J3 < —0.05 
for finite chains the transition from the fully polarized FM state with total 
spin 5* = N/2 is not directly to a singlet state with 5 = 0, rather there are 
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Figure 1: Phase diagram of the J1-J2-J3 Heisenberg chain determined by 
CCM for N -> 00 and ED for N = 24, 28, 32. Below the transition line the 
GS is the fully polarized FM state. Above the line the GS exhibits spiral 
correlations. The transition can be continuous (larger J 3 ) or discontinuous 
(smaller J3), see text. Note that the classical transition line (black dashed) 
corresponds to the result of Pimpinelli et al. [12] . 
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some intermediate states with N/2 > S > 0. To give an example, for J 3 = 
—0.08 and N = 32 there is a sequence of transitions from the FM state to a 
singlet GS via partially polarized states in the interval 0.354 < J 2 < 0.372. 

rpED rpED 

The estimation of the the tricritical point by ED is T^ U D = (J 2 qu , J 3 qu ) ~ 

(0.28, —0.01), cf. Fig. HI Note that for J 3 > J 3 qu the classical and quantum 
transition lines coincide, i.e. the relation J| = (1 — 9Js)/4 is valid also 
for the quantum model. Due to the quite large prefactor —9/4, already a 
weak 3rd neighbor coupling J3 has a noticeable effect on the critical Jf . For 
J3 < J 3 qu the transition point of the quantum model is shifted to smaller 
values of frustrating J 2 . That is to some extent surprising, since in most of 
the previous studies of models exhibiting a transition between a spiral and 
a collinear GS in the quantum model the opposite behavior has been found, 
i.e., the transition to the spiral GS is shifted to higher values of frustration, 
see, e.g., Refs. ^MM^B^M • 

To illustrate the quantum tricritical point in more detail we show in Fig. [2] 
the CCM results for the pitch angle versus J2 for various values of J 3 around 

rpLSUBA 

J 3 qu . For comparison we show also the classic pitch angle j c i. It can be 
clearly seen how the continuous behavior of the pitch angle 7 goes over into a 
discontinuous one. Interestingly, at a particular value of J2 = J| the curves 
cross each other, i.e. the pitch angle is independent of J3. For the classical 
model the crossing point is at J| = 1/2 and the corresponding pitch angle 
is lei = tc/3. For the quantum model the curves do not cross exactly in a 
point, rather they approach each other very closely at J 2 = 0.335. The pitch 
angle at that point is 7 gu = 0.327T. Furthermore the quantum pitch angle 7 9U 
approaches the limiting value limj^oo 7 gu for much smaller values of J 2 than 
the classical one. 

In Fig. [3]we present the height of the jump 7^ at the transition point in 
dependence on J 3 . For the classical model 7^ is given by Eq.(J5]) for J 3 < 
— 1/15. For the quantum model the 7t(^3) curve is characterized by two 
nearly linear regimes, one regime (near the quantum tricritical point) with 
a steep increase of 7^ and a second, almost flat one for J 3 < —0.2. This 
scenario is confirmed by the ED results for the NN spin-spin correlation 
function shown in Fig. 2] which may serve as a finite-chain analogue of the 
infinite-chain pitch angle. 

Finally, let us discuss the pitch angle which appears in the limits of large 
J2 or large | J^\. This limiting value of 7 is the maximal pitch angle and it 
is monotonously approached from below increasing the corresponding bond 
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Figure 2: The CCM results for the quantum pitch angle ^ qu in dependence 
on J 2 for various values of J3. For comparison we also show the classical 
pitch angle 7^. Note that for J3 = —0.08 (quantum and classical model) and 
for J 3 = —0.04 (quantum model only) the pitch angle jumps from zero to a 
finite value 73- at the transition point. These jumps are indicated by filled 
circles (quantum model) or a filled square (classical model). 
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J 2 or I J 3 | while fixing the other one. For J 2 — > 00 (and finite | J 3 |) the pitch 
angle approaches 7 = 7r/2. In this limit the system splits into two decoupled 
AFM chains with coupling strength J 2 . 

For J 3 — > 00 (and finite J 2 ) the pitch angle approaches 7 = 7r/3, i.e. 
only acute pitch angles appear. For J 3 — > —00 (and finite J 2 ) the pitch 
angle approaches 7 = 27r/3, i.e. by contrast to the pure Ji-J 2 model also 
obtuse pitch angles appear. As it is obvious from Fig. [3] pitch angles 7 > n/2 
appear already for quite moderate values of J3. We mention that the above 
discussion is not purely academic, since at least values J 2 > 1 might be 
realized also in real materials, e.g. in NaCu 2 2 [6]. 
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Figure 3: The height of jump of the pitch angle 7t at the transition point in 
dependence on J 3 . 
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Figure 4: The jump of the nearest-neighbor spin-spin correlation function at 
the transition from the FM to the singlet GS, cf. Fig. [TJ for finite periodic 
chains of length N=24 and N=32. 
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5 Summary 



Using the coupled-cluster method (CCM) and the Lanczos exact diagonal- 
ization technique (ED) we have studied the influence of a third-neighbor 
exchange J3 on the GS of the spin- 7; Heisenberg chain with FM NN inter- 
action J 1 and frustrating AFM NNN interaction J 2- In particular, we have 
analyzed the transition from the FM GS (present for dominating J\) to a 
singlet GS with incommensurate short-range spiral correlations. The results 
obtained by these two different approximations agree well. Moreover, the 
finite-size effects inherent in the ED study appeared to be small. 

We have found, that in case of an AFM coupling J3 the FM-spiral transi- 
tion point J% of the quantum model coincides with that of the classical model, 
and it is always continuous. However, the quantum pitch angle significantly 
deviates from the classical one. For a FM coupling J 3 quantum fluctuations 
shift the FM-spiral transition point J| to smaller values, and the transition 
becomes discontinuous. 

Acknowledgment We thank the DFG for financial support (grants 
DR269/3-1 and RI615/16-1). For the exact diagonalization J. Schulenburg's 
spinpack was used. 
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